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We report on a study of the relative rates of B meson decays into ^>(2S) and J/ip mesons using 
1.3 ftr 1 of pp collisions at i/s = 1.96 TeV recorded by the DO detector operating at the Fermilab 
Tevatron Collider. We observe the channels B° -> ip(2S)<j>, B° -> J/i>4>, B ± -> ip(2S)K ± , and 
B ± -> Jj^K^ and we measure the relative branching fractions for these channels to be 



^gt^/Sf = °- 53 ± °- 10 (stat) ± °-° 7 (syst) ± °-° 6 {b) > 
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^fe^j/fyffl = °' 63 ± 0,05 (stat) ± 0,03 (sySt) ± °-° 7 {B) ' 
where the final error corresponds to the uncertainty in the J/tp and ip(2S) branching ratio into two 
muons. 

PACS numbers: 13.25.Hw, 14.20. Mr 



B meson decays into final states containing charmo- 
nium play a crucial role in the study of CP violation and 
the precise measurement of neutral B meson mixing pa- 
rameters [T] . For CP violation in B s mixing, B s — > J/ ip<f) 
decays are being used to measure the width difference be- 
tween the mass eigenstates and the CP violating relative 
phase difference between the off-diagonal elements Ti2 
and M%2 describing the mixing of neutral B mesons [2, 3J. 
Since the current experimental results are limited by 
statistics, it is important to establish new channels like 
the decay B® — > ip(2S)(j> where these studies can be per- 
formed. 

The study of B meson decays into several charmonium 
states can also be used to constrain the long-distance pa- 
rameters associated with color octet production which 
are important for the understanding of both mixing in- 
duced and direct CP violation [4|. While these modes 
have been precisely measured in B + and B decays [5], 
an extension of these studies into the B° s system provides 
an important test of quark-hadron duality. 

In this Letter we report measurements of B meson de- 
cays into charmonium using the channels B + — *■ J/ipK + , 
B+ -> 4>(2S)K+, B° s -> J/<# and B° s -> %b{2$)4>. 
Charge conjugation is implied throughout. The J/tp and 
ip(2S) mesons are reconstructed in the dimuon channel 
and the 4> is reconstructed in the K + K~ channel. The 
study uses a data sample of pp collisions at y/s = 1.96 
TeV corresponding to an integrated luminosity of approx- 
imately 1.3 fb _1 recorded by the DO detector operating 
at the Fermilab Tevatron Collider. Similar studies have 
recently been reported by the CDF collaboration [6 . 

DO is a general purpose detector described in detail in 
Ref. [7j- Charged particles are reconstructed using a sili- 
con vertex tracker and a scintillating fiber tracker located 
inside a superconducting solenoidal coil that provides a 
magnetic field of approximately 2 T. The tracking vol- 
ume is surrounded by a LAr-U calorimeter. Muons are 
reconstructed using a spectrometer consisting of magne- 
tized iron toroids and three super-layers of proportional 
tubes and plastic trigger scintillators located outside the 
calorimeter. Only data recorded by dimuon triggers were 
used for this analysis. 

The selection requirements are determined using simu- 
lated samples for the four decay modes. The PYTHIA [5] 
Monte Carlo (MC) generator is used to model bb pro- 
duction and fragmentation, followed by evtgen [9] to 
simulate the kinematics of £?-meson decay. The detec- 
tor response is simulated using a GEANT [TO] based MC. 
Simulated events are processed through the same recon- 



struction code as used for the data. The dimuon trigger 
is modeled using a detailed simulation program incorpo- 
rating all aspects of the trigger logic. The trigger simula- 
tion is verified using a data sample collected with single 
muon triggers. Backgrounds are modeled using data in 
the mass sideband regions around the candidate B me- 
son. 

Muon candidates are required to have track segments 
reconstructed in at least two out of the three muon sys- 
tem super-layers and to be associated with a track re- 
constructed with hits in both the silicon and fiber track- 
ers. We require that the muon transverse momentum, 
Pt, is greater than 2 GeV/c. The charmonium system is 
formed by combining two oppositely charged muon can- 
didates that are associated with the same track jet [TTj 
and form a well reconstructed vertex with x 2 /DOF < 16. 
We require the dimuon pt to be greater than 4 GeV/c. 
The invariant mass of the dimuon system is required to 
be within 250 MeV/c 2 of the nominal charmonium state 
mass [5] , since the invariant mass resolution is about w 75 
MeV/c 2 . We then redetermine the muon momenta with 
a mass constraint imposed when forming the B meson 
candidate. 

All charged particles within the same track jet as the 
dimuon system are considered as kaon candidates and 
the kaon mass is assigned. The candidates are required 
to have hits in both the silicon and fiber trackers and 
have px > 0.6 GeV/c. Pairs of oppositely charged kaons 
with pr > 0.9 GeV/c are combined to form <j> candidates. 
The expected invariant mass resolution for the cj> mesons 
is 4 MeV/c 2 . Therefore the pair of kaons must form a 
well reconstructed vertex and have 1.008 < m(K + K~) < 
1.032 GeV/c 2 . 

The charmonium candidates are combined with cither 
a kaon or <fi candidate to form either a B + or a B® can- 
didate. The B meson daughter particles are required to 
form a well reconstructed vertex and have an invariant 
mass between 4.4 and 6.2 GeV/c 2 . 

Backgrounds from prompt charmonium production are 
reduced by requiring the B meson decay vertex to be dis- 
placed from the interaction point in the transverse plane 
by more than 4 times the error on the measured displace- 
ment for B + candidates and 6 times the error for B® can- 
didates. For all candidates, the error on the displacement 
measurement is required to be less than 150 /im. Com- 
binatorial backgrounds are reduced by requiring the B 
candidate momentum vector to be aligned with the posi- 
tion vector of the secondary vertex to within 26 degrees. 
Possible background contamination from kaons and pi- 
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ons misidentified as muons and other source of B decays 
which could result in a peaking background have been 
studied and found to be negligible. 

The resulting mass distributions of the B meson can- 
didates are displayed in Figs. [T]|4] The B meson yield is 
extracted from the data using a binned likelihood fit to 
the data assuming a Gaussian component for signal and 
a second-order polynominal distribution for background. 
The number of signal events is obtained by integrating 
the fit functions over the range of interest. This range is 
indicated by the two dashed vertical lines in Figs. [T]|4] for 
each channel, and covers a a region corresponding ±3a, 
where a is the expected invariant mass resolution. We 
see signals in all four channels. The results of the signal 
yields, corrected for the background contributions, are 
listed in Table [T] 
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J/ip (/> data selection. The region between two dashed vertical 
lines represents the signal window. 
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FIG. 2: fi + fj,~ K invariant mass distribution for the B — > 
ip(2S) K data selection. The region between two dashed 
vertical lines represents the signal window. 



The relative yield of B meson decays into ip(2S) and 
J/ip mesons is given by 



B{B -> ip(2S)M) 
B(B -> J/ipM) 



(1) 



N, 



i>(2S)M 



N 



J/ipM 



e J/tf>M 
£i/>(2S)M 



MM 
->■ fx+fx- 



where B is either a B + or B® meson, M is either a 
K + or meson, AT is the number of signal events re- 
turned from the fit, and e is the reconstruction effi- 
ciency determined from MC. The measured branching 
fractions B{J/ip -> = (5-93 ± 0.06) • 10" 2 and 

B(ip(2S) -> = (7-5 ± 0.8) • 10~ 3 are taken from 

Ref. [5] and combined into a ratio of branching fractions 
B{J/ip -> n + ijt) / B{ii>{2S) -> M + M~) = 8.12 ±0.89. 
The uncertainty on the ratio is given by the uncertainty 
on the single measured branching fractions assuming no 
correlations. 

For the measurement of B meson branching fractions 
the sources of systematic uncertainties are (i) the branch- 
ing fractions of the charmonium mesons to dimuons, (ii) 
systematics of the individual signal yield determinations 



and (iii) the determination of the efficiencies e 



■4>(2S)4> 



and 



e J/i/>0- I n t ne ratio many systematic uncertainties cancel, 
such as the integrated luminosity, b production and frag- 
mentation and the selection efficiencies. However, the 
polarization could be different for the B® decays. We use 
a pure CP-even state for the generated B® — > ip(2S)<fi 
and B° -► J/V>0 MC. 

The relative uncertainties that enter into the calcula- 
tion of the relative branching fractions are given in Ta- 
ble [IT] The uncertainty related to the measured charmo- 
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TABLE I: Summary of obtained event yields from the fits as 
described in the text and signal efficiencies obtained from MC 
simulations. 



Decay 




Efficiency 


Yield 


B ± ^ 
B ± ^ 


J/iPK+ 
ip(2S)K+ 


(1.07 ±0.02) • 10 _a 
(1.14 ±0.04) ■ 10" 3 


6276±97 
535±30 






(14.4 ±0.7) • 10" & 


565±26 


B° 3 ^ 


V(2S)0 


(15.2 ±0.6) • 10~ 5 


40±8 



nium resonance branching fractions enter both measure- 
ments and are the same for both. The uncertainties are 
treated as uncorrelated and give a combined uncertainty 
of 11% on each of the ratios of branching fractions. 

The relative statistical uncertainties on the efficien- 
cies e^(2S) k an d £.j/^k are 3.7% and 1.9%, respectively. 
They are combined into a single statistical uncertainty on 
the efficiency ratio assuming no correlations. To obtain 
and estimate the signal yield variation, the background 
shape and fit range for the background region are var- 
ied. This yields a variation of 3% for B ± -> tp(2S) K ± 
and 0.9% for B ± -> J /if) K ± . The ratio of branching 
fractions B(B ± -► ^(2S)K ± )/B(B ± -> J/ipK*) is then 
0.63 ± 0.05 (stat) ± 0.03 (syst) ± 0.07 (B). 

The relative uncertainties on e^su and ejM^ are 
4.5% and 5.6%, respectively. These uncertainties are 
combined into a single statistical uncertainty on the effi- 
ciency ratio assuming no correlations. As an estimate of 
the signal yield variation, the shape of the background as 
well as the invariant mass regions for the background es- 
timations are changed. This gives a variation of 7.5% for 
the B° -> ip(2S)<t> and 2.1% for the 5° -► J/</>0 signal 
yield. The B° -> J/ip^> and B° s -> ip(2S)<j> MC events 
are generated as pure CP-even decays with a B® life- 
time of 1.44 ps [12]. To account for a possible efficiency 
difference related with the different lifetime of the B®, 
the B® — > J/ip<fi MC events are weighted according to 
the combined world average lifetime 13]. The efficiency 
difference is estimated to be 8%, which is taken as an 
additional systematic uncertainty. The resulting ratio of 
branching fractions is B(B° -> ip(2S)(f)/B(B° -> J/ip<j>)) 
= 0.53 ± 0.10 (stat) ± 0.07 (syst) ± 0.06 (B). 

In summary, we have presented the observation of the 
decay B® — > ip(2S)<f) with the decay ip(2S) — > /i + /.t~ at 
DO and performed a measurement of the ratio of branch- 
ing fractions 

g(B°^V>(2S)0) 

B(B°-J/^) (> 
0.53 ± 0.10 (stat) ± 0.07 (syst) ± 0.06 (B). 

In addition, a measurement of the ratio of branching frac- 
tions 

B(B± -» J/i/>K±) { > 



TABLE II: Relative systematic uncertainties on the ratio of 
branching fractions B(B ± -> tp(2S) K ± )/B(B ± -> J/ipK*) 
and B(B° -> ip(2S)(t>)/B(B° s -> J/ip<f>). 



Source Relative Uncertainty [%] 

(gggggjg (J/y>,^(2S))0 



B(J/ip -> fj.fl) 1.7 

B(j>{2S) -> wt) 11 

Total g 11 

£ j/i/.(if±,<#>)/ e i/.(2s)(if±,(/>) 4.1 7.2 

Event yield ip(2S) channel 3.0 7.5 

Event yield J/ip channel 0.9 2.1 

CP odd-even mixture (J/i>4>) N.A. 8.0 

Total (syst) 5^2 131 



0.63 ± 0.05 (stat) ± 0.03 (syst) ± 0.07 (B) 

has been performed. These results are competitive and 
in good agreement with published measurements [HI 1X4] . 
The combination with these measurements should result 
in a significant precision improvement on the measured 
ratios of branching fractions. 

We thank the staffs at Fermilab and collaborating 
institutions, and acknowledge support from the DOE 
and NSF (USA); CEA and CNRS/IN2P3 (France); 
FASI, Rosatom and RFBR (Russia); CNPq, FAPERJ, 
FAPESP and FUNDUNESP (Brazil); DAE and DST 
(India); Colciencias (Colombia); CONACyT (Mexico); 
KRF and KOSEF (Korea); CONICET and UBACyT 
(Argentina); FOM (The Netherlands); STFC (United 
Kingdom); MSMT and GACR (Czech Republic); CRC 
Program, CFI, NSERC and WestGrid Project (Canada); 
BMBF and DFG (Germany); SFI (Ireland); The Swedish 
Research Council (Sweden); CAS and CNSF (China); 
and the Alexander von Humboldt Foundation. 



[a] Visitor from Augustana College, Sioux Falls, SD, USA. 

[b] Visitor from The University of Liverpool, Liverpool, UK. 

[c] Visitor from ICN-UNAM, Mexico City, Mexico. 

[d] Visitor from II. Physikalisches Institut, Georg- August- 
University, Gottingen, Germany. 

[e] Visitor from Helsinki Institute of Physics, Helsinki, Fin- 
land. 

[f] Visitor from Universitat Zurich, Zurich, Switzerland, 
[fl Deceased. 

[1] A. B. Carter and A. I. Sanda, Phys. Rev. D 23, 1567 
(1981); I. I. Bigi and A. I. Sanda, Nucl. Phys. B193, 85 
(1981). 

[2] T. Aaltonen et al, [CDF Collaboration], Phys. Rev. Lett. 

100, 161802 (2007). 

[3] V. M. Abazov et al, [DO Collaboration], Phys. Rev. Lett 

101, 241801 (2008). 

[4] M. Beneke, F. Maltoni, I. Z. Rothstein, Phys. Rev. D 59, 

054003 (1999). 
[5] C. Amsler et al., Physics Letters B 667, 1 (2008). 



[6] A. Abulencia et al., [CDF Collaboration], Phys. Rev. 

Lett. 96, 231801 (2006). 
[7] V. M. Abazov et al. [DO Collaboration], Nucl. Instrum. 

Meth. Res. A 565 463 (2006). 
[8] T. Sjostrand et al, Computer Phys. Commun. 135, 238 

(2001). 

[9] D. J. Lange, Nucl. Instrum. Methods Res. A 462, 152 
(2001). 

[10] R. Brun and F. Carminati, CERN Program Library Long 

Writeup W5013, 1993 (unpublished). 
[11] All charged particles in the event are clustered into jets 

using the DURHAM clustering algorithm with the pr 

cut-off parameter of 15 GeV/c. S. Catani et al, Phys. 

Lett. B 269, 432 (1991). 
[12] V. M. Abazov et al. [DO Collaboration], Phys. Rev. Lett. 

94, 042001 (2005). 
[13] Heavy Flavor Averaging Group (HFAG), arXiv:hep- 

ex/0505100] (2005). 
[14] B. Aubert et al. [BABAR Collaboration], Phys. Rev. D 

65, 032001 (2002). 



